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ABSTRACT 


A slmole, semi-empirical performance correlatlon/nredlctlon 
technique aopllcable to gaseous and liquid propellant rocket 
enqlnes Is presented. Excellent correlations were attained for 
over 100 test firings by "adjusting" the confutation of the 
gaseous mixing of an unreactive, coaxial jet using a Correlation 
Factor, F, which resulted In prediction of the experimental 
combustion efficiency for each firinq. 

Static pressure, mean velocity and turbulence Intensity In 
the developing region of non-reactive coaxial jets, tynical of 
those of coaxial Injector elements were determined. Oetalled 
profiles were obtained at twelve axial locations (extending from 
the nozzle exit for a distance of five diameters) downstream 
from a single element of the Bell Aerospace Hg/Og Id-element 
coaxial injector. These data are compared with analytical ore- 
dictions made using both eddy viscosity and turbulence kinetic 
energy mixing models and available computer codes. Comparisons 
were disappointing, demonstrating the necessity of developing 
Improved turbulence models and computational techniques before 
detailed predictions of practical coaxial free jet flows are 
attempted. 



SUMMARY 


i . 

A simple, semi-empirical performance correlatlon/predlctlon 
technique applicable to gaseous and liquid oropellant rocket 
engines is nresented. Excellent correlations were attained by 
"adjusting'' the computation of the gaseous mixing of an unreactive, 
coaxial jet using a Correlation Factor, F, which resulted in pre- 
diction of the experimental combustion efficiency for each firing. 
The technique was successfully applied to Rocketdyne, Aerojet, 

TRW, and Bell Aerospace gaseous H 2 /O 2 rocket engines utilizing 
coaxial, triolet, trislot, premix, and reverse flow Injector 
elements, and to Bell's 6000 lb-thrust Orbital Maneuvering and 
600 lb-thrust Reaction Control Engines, which utilize triplet 
and unlike doublet Injector elements, respectively, and liquid 
monomethylhydrazl ne and nitric oxide propellants. 

Test data from over 100 firings, representing a wide range 
of engine sizes and flow conditions, were compressed when the 
Correlation Factor, F, times the Chamber Length, L, divided by 
the "effective" radius, R], was plotted versus the ratio of the 
Injection velocity of the fuel divided by the Injection velocity 
of oxidizer. The resulting correlations are useful for comparing 
the effectiveness of different Injectors at the same velocity 
ratios (which are proportional to the mass flux ratios), for pre- 
dicting optimum operating conditions for a given Injector geom- 
etry, and for assessing the consistency of test data. Before 
liquid rocket and Space Shuttle Main Engine performance predictions 
can be made with confidence, additional liquid rocket data, 
covering a wider range of conditions must be .correlated. 

Static pressure, mean velocity and turbulence Intensity In 
the developing region of non-reactive coaxial jets, typical of 
those of coaxial Injector elements, are presented. Detailed 
profiles were obtained at twelve axial locations (extending from 
the nozzle exit for a distance of five diameters) downstream from 
a single element of the Bell Aerospace H 2 /O 2 19-element coaxial 
injector. Measurements of mass-flux per unit area (using a 
constant temperature anemometer); total pressure, and local 
temperature were used In the determination of local static pres- 
sure and velocity. These data show a low pressure region exists 
near the nuzzle exit. Although this pressure reduction is small 
(0.34 psl), it substantially altered the flow development. 
Comparison of results shows that velocity near the nozzle exit 
decreased initially (there was no central velocity core) as a 
result of both pressure gradients and viscous mixing. These 
data are compared with analytical predictions made using both 
eddy viscosity and turbulence kinetic energy mixing models and 
available computer codes. Comparisons were disappointing, 
demonstrating the necessity of developing improved turbulence 
models and computational techniques before detailed predictions 
of practice! coaxial free jet flows are attempted. 
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II. INTRODUCTION 

Background 

As combustors become more complex and mission considerations 
require maximum delivered performance, the need for quantitative 
prediction of performance of advanced propulsion systems becomes 
Increasingly Important. In the past, the design and optimization 
of rocket and ramjet combustors has bee* accomplished primarily 
using trial -and-error procedures. Although this approach can be 
very costly and time consuming, the extrrmely complex processes - 
for example, Injection, mixing, and combustion - that occur In a 
practical combustor, and their Interactions, have defied detailed 
analysis leaving little alternative. That is, the turbulent mixing 
reacting flow field within practical combustors, which consists of 
unsteady, three-dimensional flows, recirculation regions (at least 
In the immediate vicinity of the Injector), and significant radial 
and axial pressure gradients, cannot be predicted In detail, e.g., 
Ref. 1. Therefore, a semi -empirical approach must be "sed If 
predictions of the performance of practical combustor^ are to be 
made utilizing current limited analytical capabilities, and ex- 
isting turbulence models, e.g., Ref. 2. 

Recent success in the development of a simple empirical eddy 
viscosity model (3,4,5) suggested that it might he used In a semi- 
empirical correlation technique based on the rational assumption 
that the turbulent mixing of reactants is the rate controlling 
process in combustors employing diffusion flames. If such a cor- 
relation technique could be made to work, it would produce design 
guidance years before the extremely sophisticated computational 
techniques that are needed to characterize, in detail, the turbu- 
lent reactlnq flows of practical combustors - if indeed they can 
t be! The development of a semi -empirical correlatlon/prediction 

technique was the subject of this contract; to assist in the de- 
| velopment, various turbulence models were evaluated by comparison 

* of predictions with hot-wire anemometry data. 

Objective 

The objective of this contract was to obtain a semi -empirical 
technique fo predicting combustor performance (C* efficiency) 
for gaseous propellants. In order to accomplish this objective, 
a large number of actual rocket test data, covering a wide range 
( of test conditions and injector geometries, first had to be 

I successfully correlated. 

i 

j Program Description 

| The program consisted of applying an existing technique, 

y previously developed at Bell Aerospace Company, Ref. 6, for the 

I prediction of combustion efficiency of gaseous H 2 /Oj> rocket 
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engines utilizing the following types of injectors: coaxial, 
triplet, trislot, premix, and reverse flow. 

In order to better characterize the coaxial injector element, 
which formed the basis of the technique, detailed flow measure- 
ments were made during cold-flow tests of simulated, coaxial in- 
jectors using hot-wire instrumentation. These measurements for 
a coaxial jet in which the outer jet velocity was 700 ft/sec and 
the central jet velocity 250 ft/sec demonstrated that conventional 
modeling techniques will not apply; reductions in static pressure 
of over 0.3 psi occurred in the vicinity of the injection station, 
and that significant and unexpected radial and axial pressure 
gradients were present. 

Computations also were made using the new Spalding (BRASS) 
and the Bell Aerospace jet mixing programs, and compared with the 
data. Agreement with the actual data was not at all good. 

A binary (machine language) deck compatible with the 7090/ 
7094 Data Processing System at NASA-teRC has been prepared for 
NASA-LeRC. 

Because the correlation technique is based on actual exper- 
imental data using a semi -empiri cal approach, it appeared likely 
that liquid, and 1 iquid/gaseous-fuel ed rocket engine performance 
also could be correlated, and ultimately predicted through its 
application. The technique was successfully applied to the 
liquid-fueled Bell Aerospace Orbital Maneuvering and Reaction 
Control Engines, and the 1 iquid/gaseous-fueled Rocketdyne Space 
Shuttle Main Engine. Although beyond the scope of this contract, 
results of these analyses are included in this report because of 
their potential importance. 

III. ANALYSIS 


Basis for Correlation 

The correlation technique was based on the assumption that 
turbulent mixing alone is the rate controlling process; therefore, 
the rate of chemical reactions that occur in combustors was as- 
sumed to be very fast compared to the mixing rate. Since almost 
all practical combustors are fueled with extremely reactive pro- 
pellants, this assumption appears plausable; making it permits 
engine performance, i.e., combustion efficiency, to be computed 
without considering the finite-rate chemical kinetics, and re- 
sulted in an enormous savings in computer time. A final simpli- 
fication, which also appeared justified since only an empirical 
correlation was being sought, was that the stagnation temperature 
remains constant and equal to the initial (unreacted) propellant 
temperature before combustion, rather than using a computed 
temperature based on the extent of the mixing and reaction. 
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For the practical gaseous H 2 /O 2 injectors analyzed, excess 
hydrogen, required for cooling and performance optimization, was 
used in each of the injector elements. In such configurations, 
most of the oxygen ultimately should react with hydrogen from 
the same element; therefore, the Interaction between Injector 
elements could be neglected as a first approximation. Using this 
assumption permitted the axisymmetric mixing model and numerical 
integration procedure (discussed below) to be applied directly 
for the computation of combustor performance; that is, the calcu- 
lation of the mixing for a single coaxial injector element (or 
"equivalent" coaxial element for other injector types) could be 
assumed indicative of each of the injector elements, thereby 
permitting correlation and prediction of the performance of the 
overall combustor to be made based on the mixing of a single 
element. 

Of course, performance computed usinq a mixing model aooli- 
cable to nonreactive coflowing free jets, without pressure gradi- 
ents, could not be expected to predict actual rocket engine per- 
formance without modification. The approach selected herein was 
to empirically adjust the values of the eddy viscosity predicted 
using the mixing model of Ref. 3 (discussed below) by multiplying 
it by a Correlation Factor, F, which remained constant in each 
performance calculation, that is, the value of the turbulent 
mixing coefficient (eddy viscosity) used for the performance 
calculation became 


The value of F was obtained 
performance data until 


= Fe 


model 


(1) 


by iteration for a particular set of 


n pred w n exp 


( 2 ) 


The linearity between F and n Dre( j assisted greatly in obtaining 
the desired value of n prec j which satisfied the Inequality 


n pred 


n 


exp 


< 0.001 


(3) 


That is, the value of Op re d was within 0.1% of the value of n exp 
for all the correlations reported herein. exp 


a 



Because the Correlation Factor, F, was empirically 
determined, it includes a variety of "real -world” effects such 
as the influence of multiple-jet interaction, heat losses, the 
interaction of mixing and combustion, and local temperature and 
density variations. For example, computation of the detailed 
chemical kinetics requires instantaneous (molecular) concentra- 
tions, rather than the time-averaged concentrations. Therefore, 
even if time-averaged concentrations could be accurately com- 
puted, they could not be expected to yield appropriate reaction 
rates without some empirical modification to either the mean 
concentration or the chemical kinetic rate constant. Unfortu- 
nately, such sophisticated modeling is beyond current capabili- 
ties , e.g. , Ref. 7. 

The important feature of this approach is that it retained 
the nonlinear effects of both the partial differential equations 
(shear layer equations), and the mixing model. Using it, cor- 
relations applicable over a wider range of combustor variables 
were anticipated than when using purely empirical correlations 
alone - for example, data fitting techniques. In particular, 
better scaling was anticipated using the proposed technique than 
with previous techniques because the size of the individual in- 
jector elements and actual chamber length is used in the calcula- 
tions rather than, for example, the overall chamber diameter. 

Note, that the eddy viscosity mixing model upon which the cor- 
relation technique is based has been demonstrated to be appli- 
cable for conditions covering a reasonably wide range of geome- 
tries and velocities; data were used in the model development 
in which jet areas and velocities varied by factors of 10,000 
and 2,500, respectively, Refs. 4, 5. 

Of course, a technique in which correlation of an experi- 
mentally determined, lumped performance parameter like combustion 
efficiency is obtained assuming mixing alone to be the rate con- 
trolling process, will not predict details of the complex mixing 
and reacting flow field, such as local velocity, composition, and 
temperature. Such details only can be computed when each of the 
processes important in a combustor are included ir extremely 
sophisticated analyses that are beyond current capabilities. 
Certainly as a minimum in such analyses, the simultaneous turbu- 
lent mixing and chemical kinetics and their interactions must be 
considered, as well as the influences of multiple elements; in 
addition, recirculation effects also must be considered in such 
analyses since they provide the source of ignition once combustion 
has been initiated. 

Empirical Axisymmetric Mixing Model 

Details of the development of an empirical, mass defect, 
eddy-viscosity-type, turbulent mixing model that resulted in 
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correlation of both mean quantities and turbulence Darameters 
has been reported, Refs. 3 to 5. Because an empirical Correla- 
tion Factor, F, was to be applied for combustor correlations, 
the simpler model of Ref. 3 was selected. This model permits 
prediction of the shear stress at the velocity half width of the 
jet, which is the point at which the turbulent shear stress ap- 
proximately reaches its maximum value, and is* 


e model 


0.0180 ( 1 + 2f ) 

P 



00 


! P U-P e U e |rdr 

i 

(D-r u )exo(- 0.115F) 


(4) 


The more complex version of the model in Refs. 4 and 5 which 
predicts radial shear stress distributions as well, was not re- 
quired in the semi -empiri cal correlation technique. As a 
further simplification, no effort was made to model the core 
region of the coaxial jet separately from the transition region 
(Figure 2), even though such a distinction must be made in de- 
tailed flow analyses, Ref. 4. The relatively simple mass defect 
model of Eq. (4) was considered adequate for these correlations 
since a Correlation Factor, F, to be determined empirically, 
would compensate for the model's lack of generality, as long as 
the general functional form of the model was valid and mixing 
is indeed the rate controlling process in combustors enoloying 
diffusion flames. That is, all the complex effects which are 
present in an actual combustor and influence its performance, 
such as pressure gradients and heat losses, as well as the chem- 
ical kinetics, are lumped empirically in the single Correlation 
Factor, F. 


Combustion Efficiency Calculation 


In the correlation procedure, the combustion efficiency was 
computed by considering only the non-reactive mixinq of the H 2 
and O 2 propellants; no chemical reaction nor change in stagnation 
temperature was considered. In order to compute the combustion 
efficiency, n Dred , the assumption was made that whenever H 2 and 
O 2 were mixed Hocal ly , the HoO formed at that point could be 
computed di recti y from the limiting constituent. The total mass 
flow rate of H 2 O that resulted from this assumption was determined 
by numerical integration at each axial station; its division by 
the maximum mass flow rate of water that could be formed when all 


♦Symbols are defined in NOMENCLATURE 
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the O 2 was reacted (since H 2 always was in excess), yielded the 
predicted combustion efficiency, n pr pd . Using these assumptions, 
hpre'' increase with axial length, until ultimately it 

approaches 100%. Although this assumption obviously is not valid 
in actual combustors, since chemical equilibrium and heat losses 
present attainment of complete reaction (and hence 100% combustion) 
eve:j for very long chambers, the assumption that n D red increases 
with increasing chamber length appears reasonable tor practical 
combustors which always are reasonably short. As demonstrated 
below, these simplifying assumptions, while precluding computation 
of detailed flow conditions such as local concentration or veloc- 
ity, do not adversely affect the prediction of combustion 
eff i niency. 


The combustion efficiency was computed in the following 
manner. The local mass fraction of H 2 O, Y^O. that would have 
occurred as a result of mixing was computed as 


'h 2 o 



when Y n /Y u < 8 
°2 H 2 


(5) 


Y 


h 2 o 



when Yq /Y^ > 8 


( 6 ) 


For the cases considered, in which H 2 was always in excess, the 
maximum quant ty of water formed was proportional to the initial 
flow rate of O 2 ; that is, all the O 2 would ultimately react to 
form H o 0, so that 


\ 


^h 2 o 


) 

max 


I < G 0 > 


(7) 


Therefore, the overall combustion efficiency at any axial station 
was computed as 


pred 


2tt 

T1 V 7 Z7 

c max 


pUYn^rlr 


( 8 ) 


where roe is tin freestream boundary which defines the extent of 
the shear layer, i.e., mixing region. 

I !' v.ne computation for npred» the governing equations of 
change, which are the shear layer equations, were solved 
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numerically In von Hisses coordinates using a standard explicit I 

finite difference method, similar to that described In Ref. 8. ; 

For simplicity, the outer (external) H 2 stream was assumed to be ] 

Infinite In extent (Figure 2) which was consistent with the as- I 

sumption of excess H 2 . Of course, this meant that the velocity \ 

at the outer edge of the hydrogen stream (freestream) always . • 

remained constant, and equal to Its Initial value, which Is not 
the case In an actual coaxial injector element, as demonstrated 
in the next section of this report. Numerical integration was 
continued until the length (distance from injector to throat) of : 

the actual combustor being modeled was attained. Mean values of 
velocity, mass fractions, and static temperature were computed 
at each grid point in the flows and the value of n P red was com " 
puted from these profiles by summation at each axial station 
using Eq. (8). Computed and measured combustion efficiencies 
were then compared and an appropriate adjustment of the correlating 
parameter, F, made, and the entire numerical procedure repeated 
(starting again at the injection station) until the desired 
agreement was attained; agreement of 0.1% in n was found to be 
quite satisfactory for good correlation. Linear interpolation 
of initial results of o P red versus F reduced the number of itera- 
tions needed. The simple computer oroqram used for the correla- 
tion required only about one minute on the IBM 360/65 ner itera- 
tion; four iterations were generally sufficient for convergence. 

Additional details of the computational procedure are presented 
in APPENDIX A. 
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IV. RESULTS ANO DISCUSSION 
A. Gas Rocket Performance Predictions 

The first task of the contract required that TOO data Doints 
be selected from Refs. 9 to 13 and correlated as described in the 
preceding section. The data selected, and approved by the NASA- 
LeRC Project Manager, are presented in Tables 1 to 5. The data 
in the tables designated A were taken directly from the various 
reports (except for the central jet radius R] which was computed). 
The data in the tables designated B summarize the correlation 
obtained for each case. 

The injection velocities for the H 2 and O 2 were computed 
using the continuity equation and the perfect gas law, e.g., Uq 2 , 
the Injection velocity of the gaseous O 2 propellant was computed 
from the relation 

(60 2 )(R)(T 0 ) 

u o 2 = wjvm- 2 jr^i (9) 


A similar computation was made to calculate the injection velocity 
of the gaseous H 2 , U H . The velocity ratio, V R , was simply com- 
puted as 2 


V R 



( 10 ) 


The effective radii R] and R 2 were computed from the relations 


( 1*0 > 


TfTTET 


1/2 


«V 2 

TOT + R i 


1/2 


(ii) 


( 12 ) 


Note the thickness of the ‘splitter plate was not considered it. 
the computation of R 2 because the coaxial jet (shear layer) 
computer program cannot handle the recirculation that would be 
generated in the immediate vicinity of the splitter plate. 
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Correlation of Data 

Applying the correlation technique discussed in ANALYSIS 
using the Gaseous Hydrogen/Oxygen Rocket Performance Correlation 
and Prediction Computer Program described in APPENDIX A, re- 
sulted in the values of the Correlation Factor, F, and the Cor- 
relation Function, F ( L/R-j ) , [which is written F • L /R 1 3 tabulated 
in the last two columns of the tables designated B. Values of 
the experimental efficiency were predicted to within ±0.1% usinq 
these tabulated values. 

The experimentally determined energy release efficiency 
(or combustion efficiency) is plotted versus the velocity ratio, 
Vd, in Figs. 3A to 9A for all o* the data. The points in these 
plots show considerable scatter, which was not unexpected, since 
many different engine geometries and flow conditions were cor- 
related. The scatter of these data illustrates the magnitude 
of the correlation challenge. 

Since an empirical adjustment of the mixing was to be 
employed in the correlation technique by varying the Correlation 
Factor, F, computation of the exact experimental combustion ef- 
ficiency for each particular case always can be achieved. Of 
course, such agreement is not at all significant unless useful 
correlation of a reasonably wide range of experimental data is 
attained - that is, a useful correlation must result in a smooth, 
continuous Correlation Factor, F, over a range of test conditions 
at least for similar geometries. Therefore, the value of the 
proposed correlation technique could not be assessed until a 
considerable quantity of experimental results had been correlated. 

Figures IB to 9B are plots of the Correlation Factor, F, 
obtained for each of the injector types again plotted versus Vr. 
These figures clearly demonstrate that the chamber length, L, 
(injector to nozzle throat) is an important parameter in this 
correlation technique. Note, it is not possible to determine 
whether the velocity ratio or the mass flux ratio are the more 
basic correlating parameters since these quantities are directly 
proportional. When the pressure at the injector exit ( i . e . , P c ) 
and the temperature of each propellant is equal, the constant of 
proportionality is merely the molecular weight ratio. 

Figure 3C to 9C are plots of the final value of the turbulent 
viscosity, or turbulent mixing coefficient, e^, [defined in Eq. 
(1)3 used in the computation of the combustion efficiency; detail- 
ed discussion of the calculations are presented in APPENDIX A. 
Values of in these plots vary widely, i.e., from 0.1 to 6 x 
10"4 lbf-sec/ft2 (values of F*L/R] vary only from 3 to 17). Also, 
the e n 1 s exhibit rather irregular, erratic behavior, which clearly 
illustrates that they would be completely unsatisfactory correla- 
tion or prediction parameters. Comparison of values of e n with 
corresponding values of F (Figs. 3B to 9B), and specifically the 
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Bell Coaxial Injector results of Fig. 9C for L * 4.24 in., 
conclusively prove their unsuitability. It would be impossible 
to interpolate with confidence between any of these values in 
Fig. 9C even In the Immediate vicinity of existing points. Note 
however, that it is the magnitude of the Turbulent viscosity e n , 
(bynonymous with Eddy Viscosity and Turbulent Mixing Coefficient) 
which is the measure of the mixing that occurs in the actual 
rocket engine, not F, which merely indicates the extent of the 
" adjustment " required in emodel t e.g., see Eq. (1). 

The general similarity of the shape of the plots of F versus 
Vr (Figs. 3A to 9A) is very surprising when the great differences 
in injector geometry and operating conditions are considered. 

Only a very few points do not seem to correlate consistently, 
e.g., one of the 6-in. length Rocketdyne Coaxial injector points 
in Fig. 3B. However, the Energy Release Efficiency of this 
point, 31H, is considerably higher (2.3 and 4.2%) than adjacent 
points 25H and 33H, (Table 1A) although conditions for all three 
of them are nearly identical; therefore, the validity of this 
data point is questionable. One advantage of this correlation 
technique is it permits a simple assessment of the consistency 
of the experimental performance data. 

General Correlations 


Consideration of the Aerojet Coaxial Injector Data in Fig. 

5B suggested that F is inversely proportional to the chamber 
length, L. The correlation F«L versus Vr was applied to all the 
data and is presented in Figs. 10A and 10B. However, the values 
of F»L for the Bell Reverse Flow Injector are almost 10 times as 
high as for the other four types of Injectors (Fig. 10A) and so 
had to be plotted on a separate figure (Fig. 1 OB ) . Nevertheless, 
considerable consolidation occurred when the parameter F was 
replaced by F*L. The TRW Triplet Injector Correlation Factor, 

F, was considerably smaller than the others, as was its "effec- 
tive" central jet radius, R-j. Comparison of these results with 
those of the Coaxial Injectors suggested that F is directly pro- 
portional to R'j. Results of the F /R-j correlation of the entire 
set of data is presented in Fig. 11; again considerable consol- 
idation occurs. In fact, the curve for the TRW Triplet (R] = 
0.0256 in.) and that of the Rocketdyne Trislot (R-j * 0.105 in.) 
both of which have nearly the same length, are practically con- 
tinuous in Fig. 11; whereas, there is no indication of any cor- 
relation of these data in Fig. 10A. Results obtained with the 
Bell Reverse Flow Injector (R, = 0.713 in.) also demonstrate the 
effectiveness of the F/R ^ parameter. The maximum value of F for 
this injector is 2.2 as compared with a high of 0.32 for the 
other four types of injectors; however, values of F / R'j are only 
slightly higher for the Reverse Flow Injector than for the 
others, e.g., 3.0 vs 2.6 in." 1 , validating the usefulness of the 
F/R i parameter. 
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The final non-dimensional Correlation Function, F • L /R -j 
clotted versus Vr, fcr all the Gaseous H 9 /O 2 Rocket Engine per- 
formance Data analyzed was presented In Fig. 1. It consists of 
Correlation Factor F multiplied by l, the chamber length and 
divided by Rj, the "effective" radius of the central oxygen jet, 
and Is a combination of the F*L and F/R-j parameters. This func- 
tion compressed the performance data for all of the Injectors 
Including the Bell Reverse Flow injector; It allows direct com- 
parison of the effectiveness of various injector types at a 
given value of V R . Of course, values of F • L / R i at different 
Vr's cannot be used as a measure of the effectiveness of the 
mixing any more than can F and for the same reasons. 

Figure 1 shows that no single type of injector is best over 
the entire range of velocity ratios. Also, the functional form 
of F • L / R 1 suggests that within the range of parameters covered 
by the data, chamber length, L, is directly proportional to the 
effective radius of the central jet, R-j , for a given level of 
performance. This tentative conclusion is extremely important 
and should be verified experimentally. Nevertheless, it is 
certainly reasonable that improved mixing occurs as R-j is de- 
creased. In fact, this effect is currently utilized in chemical 
laser designs, which utilize microjet nozzles to obtain extremely 
short mixing and reaction lengths. 

The role of the Correlation Function appears to be similar 
to that of the drag coefficient for flow around bodies of various 
shapes, which predicts the shapes that will exhibit the hiqhest 
draq at a particular set of flow conditions, i.e., Reynolds 
number. The Correlation Function, F * L / R -j , permits similar com- 
parison of the relative effectiveness of various injector types 
at a particular V R because injector geometry is eliminated by 
> this simple non-dimensional parameter. Note, the Reverse Flow 

correlations were surprisingly consistent with the others in 
| Fig. 1, despite the fact that a single oxygen injector was used 

1 with R] 28 times as large as in the TRW Triplet. 

Discussion of General Correlation 

One unexpected result is the sharp break in the Correlation 
Function that occurs at Vr between 8 and 10, which very interest- 
ingly occurs very close to the overall stoichiometric ratio of 8 
at which the H ^ / O 2 completely reacts to form water. This behav- 
ior is most evident in Fig. 3B (Rocket dyne Coaxial Injector) in 
which this transition region is spanned by cases having the same 
geometry and very similar flow conditions. The most likely 
reasons for this behavior are: 

1) Failure of the mass-defect eddy viscosity mixing model 
to adequately predict details of the turbulent mixing reacting 
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flow. The model predicts mixing to be proportional to the 
quantl ty , 

It Is unlikely that a complex phenomena such as turbulent mixing 
will Increase linearly with mass defect over the entire range of 
flow conditions from Vr - 1 to 25. In addition, the density of 
the propellants at the injection station are proportional to 
their molecular weights (the molecular weight ratio is 16), and 
the velocity of the hydrogen Is always greater than that of the 
oxygen < i . e . , Vr > 1), so that a range of velocity ratios exists 
(near 16) for which very small mass defect integrals will be 
computed over a considerable length of the chamber. The singu- 
larity that occurs when this integral is zero Illustrates a 
deficiency in the model. Of course, the extent of the mixing 
that has occurred at a particular axial station influences both 
the local mean molecular weight, and hence the density, which 
both decrease with length, as well as the local mean velocity, 
which increases with length. Analysis of results plotted in 
Figs. 3C to 9C show that the values of the turbulent or eddy 
viscosity, e_, used in the actual mixing computation are generally 
near their minimum at Vr a 8, inspite of the fact that the Cor- 
relation Factors, F, are near their maximum . Since Cmode] = e n /F 
[Eq. (1)], it is clear that em 0 d e T 1s computed to be far too 
small in the range 8 < Vr < 10; nence, F must be sufficiently 
large to compensate for this defect in tie model. The plots of 
cmodel vs axial station, z, presented in APPENDIX A further il- 
lustrate this point. 

2) The hydrogen stream was assumed to be infinite in extent 
(see ANALYSIS and APPENDIX A) in the correlation calculation, 
which means that the velocity of the hydrogen stream at its outer 
edge (freestream condition) remains constant throughout the calcu- 
lations. Clearly, such a condition does not exist in a thin co- 
annular jet typical of a coaxial injector in which the velocity 
decays rapidly, as demonstrated by the experimental velocity data 
presented herein. The persistence of this high-velocity outer 
stream in the computation causes the mixing to be over predicted; 
however, the factor F, must compensate for this effect, since the 
performance of actual combustors (with a thin annular hydrogen 
jet) is being computed. 

Although the relative importance of each of these effects 
is uncertain, it is not necessary to understand them in detail, 
since accurate predictions of the engine performance (combustion 
efficiency) is obtained with the existing Correlation Function. 

The fact that practical results can be attained without detailed 
understanding of all of the physical processes involved is one 
of the major advantages of a practical semi-empirical approach 
such as the one presented herein. 
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Evaluation of Design Parameters 

The success attained with the Correlation Function when 
plotted versus the Velocity Ratio (Fig. 1), suggests that only 
those quantities explicitly Included In these parameters are 
of major Importance, and that all others are of only secondary 
Importance In the prediction of gaseous rocket engine performance. 
Two of the secondary effects discussed below are the chamber 
pressure and the film coolant level. Note, variables such as 
the injector element spacing chamber diameter, nozzle configura- 
tion, heat losses, etc., do not appear to be of prime imoortance, 
but their significance cannot be assessed directly because they 
generally are constant for each particular engine, so that there 
is no simple means for differentiating between them when comparing 
completely different engines. A test program designed soecifically 
to evaluate these secondary effects would be beneficial. Because 
of the demonstrated insensitivity of the correlation technique 
to size, this program could be conducted almost exclusively with 
smal 1 -scale engi nes . 

Those design parameters for which definite conclusions 
appear justified are summarized below: 

1. Velocity Ratio, Vr, has been demonstrated to be the 
most important parameter in the correlation (e.g.. Fig. 1). 

Note that the mass flux ratio, (pU)h 2 /(p u ) 02» 1 * * * * * * * * * * * * * 1S directly pro- 

portional to Vr and that the constant of proportionality is 

(MW^p )/ (MW o 2 ) when the static pressure and temperature of each 

propellant fs the same at the injection station. The fact that 

these ratios are important rather than the magnitude of the 

actual velocities is initially surprising. Examination of the 
mass-defect mixing model used in the correlation [Eq. (4)], and 

consideration of the Species Diffusion, Axial Momentum, and 

Continuity Equations (APPENDIX A) yields the explanation. If 

the initial values of both Ug 2 and Uh 2 are increased (or decreased) 

by the same factor, Vr, will ^remain constant ; however, the Mass 

Defect and hence e m0( j e i, [Eqs. (4) and (13)], also will be 

increased (or decreased) by this same factor (as long as the 

static temperature does not vary significantly). The Continuity 

Equation requires that the tran>vorse derivative of the radial 
velocity, and hence the radial velocity itself, increase (or 
decrease) by the same factor, because the axial and transverse 
derivatives of the axial velocity, U, both increase (or decrease) 
by the factor. Therefore, when Ug 2 and Uh? ar e varied simultan- 
eously holding Vr constant, every ^term in^the Species Diffusion 
Equation is multiplied by the factor since the turbulent mass 
mixing coefficient is computed as e«/0.7, and identical concen- 
tration profiles are obtained for all species in each case, 
which yield the same value of np re< j. However, every term in the 
Shear Layer Momentum Equation (wnen 9P/9z « 0) is multiplied by 
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the square of the factor, because the velocity throughout the 
flow field Is Increased (or decreased) by the factor as well as 
e_. Comparison of computer results demonstrated these explana- 
tions to be valid. 


Therefore, the Velocity Ratio, Vr, rather than the velocity 
of tha Individual streams Is the critical parameter influencing 


the computation of n Dre( j. The fact reasonable correlation 

of actual rocket test data was achieved with these computations, 


suggests that they adequately approximate the actual complex 
processes that occur In combustors employing gaseous fuels when 
used in conjunction with the empirical Correlation Factors. 


2. The Correlation Function, F»L/Rj, suggests that the 
ratio of the chamber length, L, to the "effective" central jet 
radius, R^, is the critical geometric parameter In the design 
of gaseous fueled rocket Injectors. This ratio Implies that 
similar performance can be expected from "scaled" combustors 
of various sizes as long as this ratio remains constant. Of 
course, such scaling can be expected to apply only over reason- 
able ranges of combustor dimensions, i.e., those within the 
range of the test data correlated herein. 

3. Chamber Pressure, P c , was shown to have little influence 
on the Correlation Factor, F, e.g., Fgs. 3B to 9B. In Fig. 9B, 

a single curve correlated the Bell 9.5-in. lenqth coaxial injector 
for P e ' s ranging from 90 to 260 psia. The overall success of the 
Correlation Function, F - L/R •) (Fig. 1) implies there is no P c 
effect, at least over the range of P c 's from 15 to above 450 psia. 
The computations are not affected by pressure because when pres- 
sure gradients are neglected, the Species Diffusion Equations, 
Shear Layer Momentum Equation, and Continuity Equation are all 
independent of pressure, as long as a mass-defect mixing model 
and the perfect gas law are used to compute emodel' Note, if 
the flow is assumed to be i soenergeti c , the momentum equation 
is not completely independent of the static temperature (as it 
is of the pressure), since static temperature depends on the 
local velocity. However, for subsonic injection velocities, 
static temperature variation is small and density variations 
are caused primarily by the molecular weight variation. 

4. The Film Coolant was shown to have little effect on the 
correlations in Figs. 1 and 3B to 9B. In fact, only in Fig. 6B 
does the Film Coolant Level appear to influence the value of F. 

As explained in ANALYSIS, all correlations were made neglecting 
the hydrogen fuel injected in the film barrier, focusing on only 
a "typical" injector element. Of course, if the combustor design 
were such that a significant portion of the film coolant were 
ingested into the combustion zone and reacted, this assumption 

no longer would be valid. Apparently, such an effect occurred 
in the Aerojet Premixed Combustor. 
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In Fig. 1 the solid squares designated Aerojet Premixed 
(No Coolant) correlate more poorly than any of th« other 100 
points. Interestingly, the Aerojet Coaxial and the Premix cases 
In which 20 to 30% of the total hydrogen was Injected as film 
coolant, all correlate reasonably well with a single straight 
line. Assuming no errors In the reported data, the most reason- 
able explanation appears to be that high levels of coolant 
Increase mixing (yielding higher values of F*l/Ri at a given Vr) 
in the Premixed Injector geometry, by causing some of the coolant 
to react, than occurs in cases In which no film coolant Is present. 
In the Premixed Injector, Ho/Op Is presumably mixed prior to 
Injection and the likelihood of the hydrogen coolant reacting 
with the oxygen present in its Immediate vicinity (near the wall) 
is obviously far greater than when each oxygen injector element 
Is completely surrounded (and hence "protected") by a.) exces' 
o f hydrogen, as is the case for most of the other tyoe . of 
injectors. Apparently, without the film coolant to enhance 
mixing, the Premixed Injector is less effective than the Coaxial 
Injector. Clearly, additional data should be correlated to 
confirm the above hypothesis. 

Application of Prediction Technique for Performance 

Optimizati on 


The ultimate application of the Performance Prediction 
Technique would be for design guidance in the selection of ortimum 
combustor geometries and flow conditions. The Correlation Func- 
tion, F*L/Ri, in Fig. 1 permits an assessment -of the Ho/Og in- 
jector configuration expected to yield the highest performance 
at the design velocity ratio. Of course, the selection of most 
efficient Injector configuration can only be made from those 
specific Injectors for which correlations already are available. 
For example, there is no means for determining whether chaogim 
the angle of impingement of a Trislot Injector element would 
improve its effectiveness - unless such an injector were test 
fired and Its combustion efficiency determined experimentally 
for at least a few velocity ratios (mixture ratios). 

Once the type of injector element to be used for a particular 
application has been selected, the fiow conditions that will yield 
the highest performance and the anticipated overa 1 1 'e'f f i ci ency 
are of interest to the designer. The Prediction Technique pre- 
sented herein may be used to select the flow conditions, e.g., 

Vr, expected to yield the highest performance. Subsequent hot 
firings would confirm the predictions and yield additional per- 
formance data, that wnen correlated, would allow refinement of 
subsequent predictions. Of course, once the velocity of each 
propellant was selected, the effective radius, R -j , of the central 
Injector could be chosen to achieve the desired mixture ratio. 
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An example of this procedure, using the Rocketdyne Coaxial 
Injector and 5-in. chamber is presented in Fig. 12, in which 
predicted Energy Release Efficiency is plotted versus Vp. In 
these predictions all input parameters including R] were held 
constant at those values indicated in the figure, except for Uo ? 
which was varied from 118 to 205 ft/sec and the Correlation 
Factor F, which was taken from the curve of Figs. 3B in which 10 
points were very consistently correlated for the Rocketdyne Co- 
axial Injector with L = 5 in. These data were selected for use 
in this illustration because the correlation curve contained a 
maximum, and the effect of this type of correlation bn the pre- 
diction was .f interest. The results in Fig. 12 indicated bcfh 
a minimum and a maximum in performance is predicted, and that 

92 to 96% - a very significant variation in 
shape of the curve is undoubtedly caused 
Correlation Factor was double valued. It 
would be extremely interesting to test fire this particular engine 
to determine whether or not both the maximum and 'minimum per- 
formance levels do indeed occur. 


n pre d varies between 
performance. The "S" 
by the C act that the 


Sensitivity of Predictions 

The Correlation Factor, F, is the only input parameter not 
read’ly determined when making performance oredictions (see 
APPENDIX A). The results presented in Fig. 1 showed that as 
anticipated in any correlation of experimental data, some scatter 
occurred about each of the correlation curves. Therefore, it 
was of interest to determine the significance of a perturbation 
in F on the predicted value of n* A specific test was made using 
the Coaxial Injection Data of Rocketdyne Case 12H by changing F 
± 4% from the value of 0.185 reported in Table IB. The ± 4% 
variation in F resulted in only a ± 0.6% variation in the pre- 
dicted value of n (ETA PRED), indicating that the variation in n 
was only about 1/6 the magnitude of the perturbation in the ore- 
dieted value of n. 

These results were typical of those obtained for ail the 
other cases, i.e., the magnitude of the slope AF/'An was very 
nearly constant for all tne cases except for the Reverse Flow 
Injectors. These results demonstrate that reasonably accurate 
pe» formance predictions can be expected when realistic estimates 
of the magnitude of F are possible, and that F does not have 
to be known precisely if estimates of n to within ± 1% are 
adequate. 
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Space Shuttle Main Engine 

Although there is little interest at present in qaseous 
H 2 /0 o rocket engines, there is considerable interest in LOX/ 
gaseous H2 engines since these propellants are used for the 
Space Shuttle Main Engine (SSME). The SSME is not only very 
much larger than any of the Reaction Control Enqines, i.e., 600 
injection elements, 8.9 in. chamber diameter and 14 in. lenqth, 
but initial designs operated at pressures as high as 3000 nsia. 
Therefore, apolying the Performance Correlation/Prediction 
Technique to this engine was a considerable challenge. Mr. 

Klaus W. Gross, NASA Marshall Space Flight Center supplied 
several sets of performance data for various enqines in various 
stages of development. Each data set was correlated by applying 
the identical technique used for the gaseous H 2 /O 2 engines, but 
the injection velocity of the oxygen was computed assuming the 
oxygen to be gaseous rather than liquid (see LIQUID ROCKET 
PERFORMANCE PREDICTIONS). 

Detailed results will not be Dresented because extremely 
limited test data were available at the time the correlations 
were made. However, values of F, F*L, F/R-j, and F*L/R, were 
all within the range of the data plotted in Figs. 1, 10, and 11. 
Obviously, definite conclusions cannot be drawn from these 
preliminary data; however, it is clear that SSME data can be 
successfully correlated, and meaningful predictions made, once 
a data base has been established. The fact that P c and the 
number of injection elements both were increased nearly an order 
of magnitude from previous maximum values made little difference 
in the magnitude of the correlation parameters. 


B. Liquid Rocket Performance Predictions 

The successful application of the semi -empi rical correla- 
tion/prediction technique to gaseous H 2 / 0 2 rockets suggested 
that it might be applied as well to liquid rockets. Of course, 
in liquid/liquid rocket engines, droplet formation, injection, 
penetration, spreading, mixing, burning, and vaporization occur 
as well as gaseous mixing. Therefore, the magnitude of the 
"equivalent mixing" of a cold-flow, coaxiai gas/gas injector 
element would be expected to be considerably smaller when applied 
to liquid/liquid injectors than when applied to gas/gas injectors. 
However, the correlation/prediction technique, because of its 
semi-empirical nature, still should be applicable. 

To demonstrate this capability, Bell Aerospace Reaction 
Control Engine (RCE), Ref. 14, and Orbital Maneuvering Engine 
(OME), Ref. 15, test data were correlated; in each case, the 
fuel was monomethyl -hydrazine (MMH) and the oxidizer was nitric 
oxide ( N 2 O 4 ) . However, two very different types of injector 
configurations were used in these engines; the RCE elements were 
unlike doublets, while the OME elements were triplets (two MMH 
jets impinging at 32 c on a central N^O^ jet). 

In order to demonstrate the applicability of the correla- 
tion/prediction technique to liquid propellants as simply as 
possible (without any modification of the computer program), 
the calculation of propellant "injection velocities" were made 
assuming that the propellants were completely vaporized, so 
that the injection velocities UmMH and Un 2 04 could be calculated 
combining the Continuity Equation and Perfect Gas Law as in Eq. 
(9), but using the appropriate molecular weights. Therefore, 
the velocities computed in this manner are referred to as "pseudo" 
velocities, since they are purely fictitious. Of course, the 
actual injection velocity of liquid propellants depends on the 
density of the liquid and not their molecular weights; however, 
computation of "pseudo" velocities was the simplest way to apply 
the technique to liquid rocket data. Since the demonstration 
calculations were made with the same computer program used for 
the H 2 /O? propellants, i ncorrect molecular weights (i.e., 2.016 
and 32.0) and thermodynamic properties were associated with MMH 
and N 2 O 4 in the n P red computations. Nevertheless, because these 
computations all were made in a consistent manner, the Correla- 
tion Factor, F, "adjusted" the mixing [Eq. (1)] to the appropriate 
level for each experimental firing and excellent correlation of 
all 21 cases was obtained. 

The experimental combustion efficiency data, n P red> is 
plotted versus the Mixture Ratio, N 2 0 4 /MMH, in Fig. 13a. The 
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Mixture Ratio was plotted from right to left so that trends would 
correspond roughly to the previous plots in which Vr was used as 
the independent variable. Clearly, attempting to correlate such 
widely scattered data was a considerable challenge for the cor- 
relation/prediction technique. Resulting correlations are ore- 
sented in Fig. 13b in which the Correlation Factor, F, is also 
plotted versus the Mixture Ratio, N 2 O 4 /MMH; these correlations 
were unexpectedly good. All the points for each type of injector 
fall on a smooth continuous curve, suggesting that meaninqful 
predictions may be made for liquid propellants as well as for 
gaseous propellants. Of course, each of the injectors has a very 
different correlation characteristic because two very different 
configurations were used. Naturally, a great deal of additional 
liquid rocket performance data must be correlated for various 
types of injectors and ranges of flow conditions before predictions 
can be made except within the rather limited range of the Bell 
OME and RCE data. Nevertheless, these results suggest that the 
relative merits of the unlike doublet and triplet be carefully 
evaluated in subsequent engine development programs. 

The final values of the Turbulent Viscosity, e n , comouted 
at x = L for the Liquid/Liquid Rocket data are presented in Fiq. 

14. These results confirm the conclusions concerning Figs. 3 C 
to 9C , i.e., the mixing coefficient itself cannot be used for 
correlation or prediction of rocket performance data. Of course, 
some of the variation in Fig. 14 is probably caused by the fact 
that "pseudo" velocities rather than actual liquid velocities 
were used. For this same reason, it is not Dossible to compare 
the magnitude of the F's for gasecus and liquid injectors and 
draw any definite conclusions. 
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C. EXPERIMENTAL COAXIAL JET FLOW CHARACTERIZATION 

The Gaseous Rocket Performance Correlation/Prediction Technique presented herein is based 
on the ability to predict, at least approximately, the non-reactive turbulent mixing that occurs in a 
coaxial jet. Therefore, it was of considerable interest to characterize such a flow in detail, and to com- 
pare the results with predictions made using available computational techniques. In this task, the Bell 
19-element Coaxial Injector, which had been previously tested during hot firings was selected to insure 
that this investigation would be relevant to practical rocket combustors. For simplicity, all elements 
except for the central element, were masked so that the complex 3-D interaction between elements 
need not be considered. A photograph of this injector is presented in Fig. 15. A schematic diagram 
showing the dimensions of the element is shown in Fig. 16. Note, that the splitter plate separating the 
central and annular jets was recessed 0.03 in. from the Sector face. 

The motivation for this investigation was twofold: First, because coaxial elements form tnc 
basic injc tion mechanism for many combustors, it was of interest to better cha» .icterize the initial flow 
region of such elements. Second, since computations for reactive systems in the past have often been 
based on the ability to predict the non-reactive turbulent mixing, it was of interest to characterize such 
a flow in detail and compare the results with predictions made using available computational techniques. 

The measurements presented herein are applicable to the developing region of a coaxial jet in 
which the injection velocity of the central jet is significantly less than that of the external annular jet. 
That is, an effort was made to characterize the classical core and transition regions (Fig. 2) of a practi- 
cal Coaxial Injector Element, in which the velocity ratio, V^, was 2.8, a value typical of the Coaxial 
Injectors for which the Correlation/Prediction Technique was applied, e.g.. Fig. 1 and Tables 1 , 2, and 
5. 

Several previous experimental investigations have been conducted to characterize the flow 
field of coaxial jets. Determinations of the mean velocity field obtained using pitot probes are dis- 
cussed in Refs. 16 to 18. Two more recent investigations have been conducted using hot wire anemo- 
meters in which both the mean and turbulence structure of the flow were determined. Refs. 19 and 
20. One assumption, which seems to be common to these investigations was that static pressure gradi- 
ents in the developing region of the flow are small and could be neglected. However, as demonstrated 
herein, when the velocity ratios. U Q /Uj, are greater than one, and area ratios, Aj_/Aj, greater than zero 
(typical of real combustors), neglecting differences in static pressure in the vicinity of the injector may 
not be a valid assumption. When using total and static pressure data to calculate mean velocity for in- 
compressible flows, even extremely small changes in static pressure can cause large errors in the velocity 
determined. These small changes in static pressure do not, however, cause large errors in velocity 
measured using a hot wire. Using hot wire anemometry, both Champagne and Wygnanski, (19) and 
Durao and Whitelaw, (20), investigated incompressible coaxial flows where the velocity ratios were 
greater than one. The measurements of Champagne and Wygnanski were concerned with the developing 
region of coaxial jets where the initial flows were nearly fully potential and separated by a very thin 
splitter plate. In contrast, the measurements of Durao and Whitelaw were for fully developed flows 
with a significant splitter plate thickness. 

The present study was concerned with the flow from a coaxial injector element from an actual 
combustor As with most practical flow hardware, the nozzle and annulus flows were neither fully 
potential nor fully developed and the separation distance (splitter plate) was relatively large. 
Measurements were made for a velocity ratio greater than one and in contrast to previous investiga- 
tions the velocities were in the compressible range. As noted earlier, these higher velocities are typical 
of real combustors. 
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Experimental P rocedure 

In this experimental investigation, a Bell 19-element coaxial injector, which had been pre- 
viously tested during hot firings, was selected and all elements except the central one were masked, 

(1 3). A photograph of this injector was presented in Fig. 1 5 and geometrical shape and physical size 
are shown in Fig. 16. Note, that the splitter plate separating the central jet and annulus was recessed 
0.03 in. from the injector face. The constant area inlet length was about 6 diameters for the central 
nozzle and about 7 nozzle thicknesses for the annulus. Therefore, fully developed turbulent flow 
conditions were not expected at the injector exit. The area ratio between the outer and inner 
nozzles, A 0 /Aj, was 1.48 and between the splitter plate and central nozzles, A l /Aj. was 1.02. The 
ratio of the annulus velocity (0 o = 700 ft/sec) to that of the central jet (U ; = 250 ft/sec) was equal 
to 2.8 for all tests. These velocities and dimensions correspond to Reynolds numbers of 3.8 x 
10 s based on diameter for the inner jet and 2.0 x 10 6 based on annulus height for the outer jet and 
are approximately an order of magnitude larger than those investigated in Ref. 19 and 20. 

In these experiments, clean, dry compressed air was supplied at constant pressure and tem- 
perature to the two plenum chambers. Plenum pressure was measured on a mercury manometer 
for the high velocity outer jet and on a precision pressure indicator (T.1.-145) for the low velocity 
inner jet. Both plenum and room temperature were maintained at 70± 1°F during all tests. 

Measurements of mean velocity and turbulence intensity were made using a linearized con- 
stant temperature anemometer (TSI-1050). The hot film probe (TS1- 1270-10) was calibrated and 
linearized to 2% of reading over the range 100 to 700 ft/sec. Below 100 ft/sec the error increased 
becoming about 4 r /< at 50 ft/sec. To assure calibration accuracy, the zero and maximum velocity 
points were checked before measurements were made at each axial station. The mean value of the 
linearized signal was determined using a true integrator (DISA-52B30) and displayed on a digital 
voltmeter (Fluke-8200A). The rms signal was obtained from a true rms meter (TS1-1060). 

Total pressure measurements were made using a pitot probe fabricated from a 0.020 in. 
outside diameter by 0.004 in. wall stainless steel tube with the tip flattened to 0.006 in. inside and 
ground to 0.010 in. outside height. Pitot pressure was sensed using a precision pressure indicator 
(Tl-145) which was calibrated and shown to have inaccuracies less than ±0.001 psid at the 3 sigma 
level of confidence. Atmospheric pressure was measured using a similar instrument which had in- 
accuracies less than ±0.004 psiu at the 3 sigma level. 

Surveys of the flow field were made by mounting the coaxial injector on a translating table 
(mill bed) and moving the injector relative to a fixed probe position. By driving the gears of the 
translating table in only one direction, the displacement of the test nozzle was known accurately 
to ±0.001 in. Initial location of the probes relative to the coaxial injector was determined using a 
theodolite. 

Detailed velocity and pressure measurements were made along the jet centerline beginning 
0.010 in. downstream of the nozzle exit and continuing for 2.50 in. Axial stations at which velocity 
profiles were obtained are 0.010. 0.030, 0.075.0.15,0.25,0.50.0.75, 1.00, 1.25. 1.50, 2.00. and 
2.50 inches. Pressure profiles were obtained at three axial stations: 0.010. 0.10. and 0.25 inches. 

Qualitative Description of Coaxial Jet Mixing 

It is helpful in understanding these results to form a picture of the developing flow field. 
Some of the flow parameters which must be considered are Reynolds numbers, velocity and area 
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ratios, and state of the How. i.e., fully viscous or potential like flow. In forming a picture of the 
flow, first consider conditions where the Reynolds numbers and velocity ratio are high (Re 0 > 

Rej > 10 4 , U D /Uj > 1 ) and where the area of the splitter plate to that of the central nozzle is ex- 
tremely small, i.e.. where the thickness of the plate separating the inner from the outer flow can be 
neglected, A L /Aj * 0. For these conditions and completely viscous flow at the injector exit, veloc- 
ity at the nozzle centerline must increase due to viscous mixing between the high velocity outer jet 
and the low velocity inner jet. However, for a condition of potential like flow -t the injector exit, 
velocity along the centerline will remain constant until viscous mixing at the boundary between the 
‘wo jets has had time to reach the centerline. In other words, a potential core will initially exist 
the length of which is dependent on both velocity and area ratios as well as the initial thickness of 
the boundary layers at the nozzle exit. An example of this type flow is discussed in the paper by 
Champagne and Wygnanski (19). 

Now consider for a moment the same flow conditions as above except let the splitter plate 
have a finite thickness, A L /Aj =% 1 . The flow must now expand to fill the region immediately down- 
stream of this plate. Therefore, a recirculation region will develop which contains two vortices of 
opposite direction of rotation, see Fig. 17. This recirculation region has approximately the same 
thickness as the splitter plate at the nozzle exit and tapers to zero at some downstream point. For 
the condition of a thick splitter plate, as opposed to a negligible one, the momentum transfer is 
initially from both jets into the recirculation region instead of from the high speed outer jet directly 
to the low speed inner one. Thus, for completely viscous flow at the nozzle exit the centerline 
velocity will decrease. Downstream where the recirculation zone ends, the momentum exchange is 
again from the outer to the inner jet and the centerline velocity will now begin to increase. An 
example of fully developed jet mixing with a finite splitter plate is discussed in the paper by Durao 
and Whitelaw (20). 

The velocity decay in the early developing region of the flow, the region contained within 
the recirculation zone, is caused by at least two separate mechanisms: viscous mixing and pressure 
gradients. To help in understanding the separate effects each of these mechanisms has on the flow 
development, consider again the flow condition shown in Fig. 17, but in contrast to the earlier ex- 
ample consider the flow at the nozzle exit to be potential with only a thin viscous region at its 
boundary. If the static pressure were constant everywhere downstream of the nozzle, as in a free 
jet. then viscous mixing would cause a decrease in velocity near the jet edge but a potential core 
would exist near the center, i.e., the centerline velocity would remain constant for a short distance 
downstream. The static pressure is. however, not constant throughout the flow. For conditions 
where the splitter plate is not negligible, the recirculation zone forms a pseudo-diffuser, causing a 
low pressure region to exist near the nozzle exit. Within this region lines of constant mass diverge 
and even for potential like flow velocity decreases almost uniformly across the jet. Viscous effects, 
which are always present, cause a slightly more rapid decrease in velocity at the jet boundary. For 
potential flow, where the viscous mixing region at the jet boundary has not spread to the centerline 
by the time the flow reaches the end of the recirculation zone, the velocity at the centerline will 
after leaving the zone, remain constant until viscous effects reach the center. 

Experimental Results 

Due to physical limitations. Hows in practical injectors are neither 
completely potential nor fully developed at tin* nozzle exit. Experimentally determined initial 
velocity profiles for one element of the Bell Id-element injector are shown in Fig. 18. Shown in this 
same figure are power law curves representing a state of fully viscous initial flow. Empirical fits, 
using experimental data, show that a I /7th power law gives the best agreement for Reynolds num- 
bers near 3.8 x 10 s and that a 1/1 0th power works best for Reynolds numbers near 2.0 x 10 < ’.(21). 
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A comparison of these curves with the measurements indicates that the flow was not fully developed. 

Note that a small potential flow region exists near the center of each jet. 

An overlay of velocity profiles measured at twelve axial locations is shown in Fig. 19a. Before dis- 
cussing these data, several comments should be made regarding the hot film used and its sensitivity. 

First, the power dissipated by the air passing over the heated film depends on both the mass flow 
rate past the sensor and the temperature difference between the sensor and fluid. Where this temp- 
erature difference is held constant, as it was during these measurements, the hot film is primarily 
sensitive to the fluid mass flow per unit area pU. For measurements made at the same flow conditions 
for which a sensor is calibrated, velocity can be determined directly; however, for measurements made 
at conditions other than those for which it was calibrated, a correction for density must be made. 

The actual velocity is thus equal to the indicated velocity times the density during calibration divided 
by the density of the flow being measured, U = Vp c /p. Second, cylindrical hot films are sensitive to 
flow velocity normal to their axis and almost insensitive to any velocity components along their axis. 

These films are, however, not sensitive to the direction of flow normal to their axis, i.e., the magni- 
tude of the velocity component normal to the axis is measured and not its direction. Thus, measure- 
ments in regions of flow recirculation will be incorrect. 

The profiles presented in Fig. 19a show a rapid velocity decrease in the initial region of the 
central jet owing to both viscous and pressure effects. Since the initial velocity at the nozzle center- 
line can be considered potential, see Fig. 18, the decrease in velocity between stations z * 0.010 in. 
and 0.075 in. must be to a large extent due to pressure effects. Evidence of strong viscous effects 
can also be seen by noting the much more rapid decrease in velocity at the jet boundary when 
going from z = 0.01 to 0.15 in. This figure indicates that for z < 0.25 in. the momentum exchange 
is from both the inner and outer jets into the recirculation zone behind the splitter plate. Some- 
where between z - 0.25 and z * 0.50 in. the recirculation zone ends and the momentum transfer is 
from the high velocity annulus flow directly to the lower velocity inner jet. Also, illustrated by this 
figure is the tendency for the higher speed annulus flow to collapse around the central jet. A com- 
parison of the peak velocity at z = 0.1 5. 0.25. and 0.50 in. shows these peaks to be moving toward 
the nozzle centerline. Figure 19b shows this trend to continue until the flow becomes jet-like some- 
where between z = 1 .5 and z = 2.0 in. It is also evident from this figure that the velocity near the 
outer edge of the outer flow, r *= 0.22 in., remains almost constant until the flow becomes jet-like, 
at which time the velocity across the entire flow field decays with distance. This is seen by noting 
that the profiles measured at all stations except z = 2.5 in., cross at approximately the same point. 

In order to quantify the importance of pressure gradients in the initial regions of the flow, 
a measure of static pressure was needed. Because direct measurements of static pressure are diffi- 
cult, if not impossible to make near the exit of such a small injector, an indirect measure was made. 

In contrast to static pressure, total pressure can be accurately measured for the flow conditions 
considered here. Measurements of total pressure, P 0 . in addition to mass flux per unit area. pO, 
and total temperature, T 0 , are sufficient to evaluate local static pressure. Static pressure at a 
point can for an ideal gas be written in terms of the local total pressure and Mach number as 

_L_ 

P„ ri , ■»-' 

~= <1+ — M") (14) 

where Mach number is defined by 
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Mach number can be rewritten in terms of velocity measured using the hot film, U - p L ,V/p, and 
static pressure as 


M 2 =~ f “~) 2 

tRT 0 \ P / 2 


(15) 


where the How temperature is equal to the temperature during calibration. Equations 14 and 15 can 
now be simultaneously solved for static pressure. 

The effect of experimental error when evaluating static pressure can most easily be seen 
by expanding Eq. 14 for small Mach number. 


P o _ ?M 2 tM^ + 

P + 2 + 1T + 


(16) 


By substituting Eq. 15 for m 2 one gets a quadratic equation for static pressure, which when solved 
becomes 


-f j p 0 V p ° ' i 


2P c -2 T*1 - tM 4 

V (l+~ M 2 )(1 + -j- 


(17) 


During this investigation, the central jet Mach number was always less than 0.23, therefore, the 
terms ( 7 -I ) M 2/2 and 7 M 4/8 were neglected and the above equation differentiated to give a simple 
expression relating static pressure error to errors in measurements of mean velocity and total 
pressure. 
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Evaluation of Eq. 1 8 for errors of 2% in V and 0.25% in P Q show the expected error in static pres- 
sure to be less than 0.4%. 

Using the hot film and pitot tube data and Eq. 14 and 15, static pressure was determined. 
The data plotted in Fig. 20 shows static pressure to be approximately 2.5% below ambient at the 
nozzle exit (14.30 psia compared to 14.64 psia), to increase rapidly and become equal to ambient 
just beyond z = 0.25 in. (approximately one inner nozzle diameter). Shown in Fig. 21 are three 
radial pressure profiles The profiles measured at z * 0 01 . 0 10 and 0.25 in show static pressure 
across the central jet to be practically constant. There does appear to be a very slight decrease in 
pressure near the jet boundary for / = 0.01 in. These figures show axial pressure gradients to be 
much stronger than the radial ones for the flow contained within the central jet. However, outside 
this region, in the recirculation z^nc, radial pressure gradients must be very steep because it is 
across this very narrow region that the pressure increases to become ambient. 
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Disc ussi on 

Figure 22 has been prepared to show that the static pressure measured at the nozzle exit is 
of the magnitude one should expect. Shown in this figure is a plot of inner nozzle exit pressure as 
a function of annulus plenum pressure. With no flow in the central jet, increasing velocity in the 
annulus causes a nearly linear decrease in pressure at the nozzle exit. Also shdwn in this figure is 
one static pressure point measured for a inner jet velocity of 250 ft/sec and an annulus velocity of 
700 ft/sec (plenum pressure 10 in. Hg). As expected, the pressure was lowest with no flow in 
the nozzle and increased with increasing velocity. 

Going back to the idea that flow downstream of the splitter plate forms a pseudo-diffuser, 
one can be integration of the data given in Figs. 19 and 20 establish contours (lines) within which mass 
is constant and equal to the total flow through the central nozzle. Doing this, one finds that the 
flow expands approximately 60% of the distance across the splitter plate at the location at which 
pressure gradients disappear. 


Effects of neglecting pressure gradients when evaluating velocity from hot film and pitot 
tube data can be seen in Fig. 23. Three separate curves of centerline velocity are shown. The velocity 
curve indicated by open circles is evaluated from hot film data directly while the curve indicated by 
solid circles is the same data corrected for changes in density (true velocity). Comparison of these 
curves shows the maximum error in velocity due to neglecting density variations to be less than 3%. 
I he third curve, indicated by diamonds, shows the velocity calculated directly from the total pres- 
sure data assuming static pressure everywhere equal to ambient (as is often done in reduction of 
pitot tube data for coaxial jets). Determination of velocity in this way can result in tremendous 
errors, e.g., the 2.5% difference in static pressure results in a difference of approximately 40% in 
calculated velocity. This large error results from taking the difference between pressure ratio and 
one in contrast to hot wire errors which results from pressure ratio only. 

Figure 23 shows not only a rapid decrease in centerline velocity due to both viscous and pres- 
sure effects (z < 0.35 in.), it also indicates a short region beyond the recirculation zone where the 
velocity remains constant (0.35 in. < z < 0.55 in.). The reason for this, as has already been ex- 
plained. is that the growing viscous region has not spread to the centerline before the flow leaves 
the recirculation zone. The sharp increase in velocity for z > 0.55 in. shows that viscous effects 
have finally reached the centerline. Beyond z 0.35 in. and until the flow becomes jet like at 
z =*= 2.0 in., the momentum exchange is directly from the high velocity annulus flow to the lower 
velocity nozzle flow. There is. of course, during the entire flow development, a viscous region 
growing at the outer edge of the annulus. Where the ratio of annulus area to central nozzle area. 
Ao/Aj, is less than one. this viscous mixing region will have a large effect on the flow development. 

For Ao/Aj > I . which is the condition considered in this paper, the effect of this viscous region is 
to cause the flow to become jet-like more rapidly (19). 


To demonstrate the accuracy of the velocity measurements, momentum was calculated at 
each axial measuring station. The percent deviation in momentum from the mean ot all measuring 
stations is shown as a function of axial displacement in Fig. 24. The maximum error over the en- 
tire measuring region is less than 2','. Two factors contributing to this small error are the non- 
lineaiity of the hot film response and the finite legion over which the momentum integral could be 
evaluated (the region over which measurements were made). Using this same data, the rate of mass 
entrainment was determined. It is apparent from Fig. 25 that the rate of entrainment remains al- 
most constant with downstream distance. Comparing this rate with that ol a tree |et. e.g.. see 
Spalding (22), the coaxial injector corresponds to a free jet with a nozzle diameter of 0.27 in. This 
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is smaller than the annulus diameter (0.44 in.) and larger than the central nozzle diameter (0.24 in.), 
thus, the measured entrainment rate is within the expected range. 

Ip addition to mean velocity, axial turbulence intensity was determined at each measuring 
station. Figure 26 shows an overlay of turbulence intensity profiles measured at seven axial stations. 
During the early development, three intensity peaks exist, two corresponding to the mixing regions 
between the recirculation zone and the jets and one corresponding to the mixing region at the outer 
edge of the annulus. The two peaks downstream of the splitter plate rapidly merge and become one 
at the end of the recirculation zone. Note that beyond this zone the inner peak decays and moves 
towards the centerline, becoming jet-like between z * 1.5 and 2.0 in. In contrast the outer peak re- 
mains at a fixed position, z * 0.22 in., while decreasing in magnitude (as for a free jet). The change in 
magnitude of turbulence intensity along the jet centerline can best be seen from Figure 27. Turbulence 
intensity increases almost linear with distance from the nozzle exit to the end of the region of re- 
circulation. z < 0.25 in. For a short distance outside this region. 0.25 < z < 0.45 in., the intensity 
remains fairly constant then increases similar to the mean velocity. These regions of similarity, 
between the centerline velocity and turbulence intensity, do not coincide \.i»n one another as can be 
seen from Figure 28. Shown in this figure are three curves - one for turbulence intensity, one for 
mean velocity and a third for the ratio of turbulence intensity and mean velocity (relative turbulence 
intensity). Each of these curves show regions of rapid change in the flow properties with axial dis- 
placement. Notice, however, that these regions are different in length as well as being slightly dis- 
placed from curve to curve. This is not surprising since turbulence and mean time scales are different. 

Figure 29 has been prepared to summarize in a concise way some of the key features of 
coaxial jet flows where the splitter plate has a finite thickness and where the annulus velocity is 
greater than the centtal jet velocity. Shown in this figure are measured velocity profiles plotted as a 
function of distance from the injector exit. The key features illustrated by this two-dimensional 
representation of the flow field are: first, the recirculation zone downstream of the splitter plate, 
second the rapid decrease in velocity of the central jet in this recirculation region owing to both viscous 
and pressure effects, and third, the eventual collapsing of the annulus flow around the central jet. 
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0. Analytical Coaxial Jet Flow Characterization 

The ability to make meaningful calculations of the 
performance of rockets and ramjets has been an Important goal of 
Propulsion Engineers for many years. Within the past five years, 
Chemical Lasers have reinforced this Interest, because It has 
become Increasingly apparent that adequate laser design requires 
the ability to predict In detail both laminar and turbulent flows 
(including recirculation regions), together with the chemical 
kinetics, and their complex Interactions. At present, adequate 
design calculations cannot be made using existing comDutatl ona 1 
and modeling techniques, e.g., Ref. 2. 

Historically, Investigators have been concerned with very 
simple flows such as single-component, non-reacting boundary 
layers, and axisymmetrlc free jets. Of course, when the com- 
plexities of turbulence must be included in analyses, the nre- 
dlction of even such relatively simple flows becomes a formidable 
computational task. Thes? complexities are compounded when 
boundary conditions applicable to actual hardware must be 
considered . 

In the Analytical Coaxial Jet Flow Characterization nhase 
of the contract, computations were made uslnq state-of-the art 
comoutational and modeling techniques. A simole coaxial qeometry 
was being comouted In which no significant density variations 
occurred. There are two different approaches qenerall.v followed 
for such analyses of turbulent flows: a) an eddy viscosity 
turbulence model, e.g., Refs. 5, 23, 24, and b) a turbulence 
kinetic energy turbulence model, e.g., Refs. 25 to 27 and 29. 
Because operational computer programs were available at Bell 
Aerospace which used the eddy viscosity model of Ref. 5 and the 
turbulence kinetic energy approach of Spalding et al, e.a., 

Ref. 29, these specific models were used for the prediction of 
the detailed axial velocity distributions for direct comparison 
with the experimental data previously presented (Section C). 

Note that the eddy viscosity model of Ref. 5 had been specifically 
developed for predicting axisymmetric coflowing streams and, 
therefore, was expected to make predictions as good as any such 
model . 

The computational technique used for the numerical integra- 
tion of the continuity, momentum, and energy equations was 
similar to the method used by Zeiberg and Bleich, Ref. 8, i.e., 
the equations were simplified using boundary layer approximations, 
the radial velocity was eliminated through use uf the von Mlses 
coordinate transformation, and the resulting equations solved by 
employing an explicit finite difference method. The turbulencp 
kinetic enerqy (TKE) nrogram and TKE model had been develooed by 
Spalding and coworkers for predicting flows bounded by walls. 

Refs. 28 and 29. In applyinq this oroqram to Predict coaxial 
free jet flows, the walls were moved from the coaxial element as 
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far as possible to minimize their effect on the computations. 

The predictive technique used, was to transform the separate 
conservation equations (mass, momentum and energy) Into a 
"standard" set of nonlinear, elliptic, partial differential 
equations, Ref. 28. The resulting coupled equation set was 
solved by a finite-difference procedure, using the Gauss-Seldel 
Interactive technique. 

Comparison of the predictions obtained with each of these 
mixing models and the experimental data (Section C) are presented 
In Figs. 30 to 41. For the eddy viscosity model, predictions 
were available at the Identical axial stations for which data 
were taken; however, because of the aspect ratio constraints of 
the grid In the "Spalding" analysis. Ref. 28, Identical stations 
were not available In all cases. In such cases, the axial sta- 
tion closest to that desired was presented as Indicated In the 
figures. Clearly, the agreement obtained with the TKE model and 
the data would not have been improved by a minor adjustment of 
the grid spacing. 

Both type of predictions were made assuming initial slug 
(step) profiles at z*0 using bulk mean velocities, as discussed 
In Ref. 4. These comparisons clearly demonstrate that neither 
type of mixing model Is applicable to the coaxial jet flow of a 
practical coaxial injector. (Features of this flow were discussed 
in some detail in Section C.) In this particular Instance the 
eddy viscosity mixing model does a far better job than the TKE 
model, although the mixing rate Is underpredicted somewhat by 
the eddy viscosity model but drastically overpredicted by the 
TKE model. Note that neither of the analyses account for the 
transverse pressure gradients that occur In the actual flow which 
are In part responsible for their failure to predict the detailed 
behavior. 

The eddy viscosity model made no adjustment for the core 
region, as was found necessary in predicting the data presented 
In Ref. 4. Had a decreased value of the eddy viscosity been used 
In the core, a decrease in the predicted rate of mixing would 
have occurred, and agreement with the experimental data would 
have been even poorer than presented in Figs. 30 to 41. 

Two important features of the eddy viscosity predictions 
are: 1) The decrease In the central jet velocity at z < 0.5 in. 

Is not predicted. Since only momentum transfer between the low- 
velocity central jet and high-velocity outer jet is considered, 
no velocities less than that of the original central jet can be 
computed, and 2) Mixing in the vicinity of the centerline Is not 
sufficiently fast for the velocity maximum to reach the center- 
line between z » 1.5 and 2.0 In. as occurs In the data. In fact, 
the predictions do not even exhibit a local maximum at the 
centerline at the 2.5 In. station. The eddy viscosity model of 
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Ref. 5 had to be multiplied by a factor of 1.7 before the maximum 
velocity occurred at the centerline at z « 2.5 In. However, In 
these predictions mixing at most radial positions were much too 
fast, and jet spreading was predicted to be considerably too 
great . 

Another eddy viscosity model that was operational at Bell 
was that of Ref. 30. Therefore, it was of Interest to compare 
predictions of this model with that of the model of Ref. 5, that 
had been under development for several years, e.p., Ref. 3. The 
predictions made with the model of Ref. 30 were considerably 
poorer than those using the original model; mixing -ates were 
still further underpredicted; therefore, detailed results are 
not presented. At the z * 2.5 in. station, the velocity profile 
still exhibited a local maximum (511 ft/sec) at a radial position 
of 0.12 In. while the velocity at the centerline was 163 ft/sec 
lower ( 348 ft/sec). Therefore, none of the profiles shapes were 
as valid as those obtained using the model of Re i . 5. Clearly, 
the model of Ref. 30 is not an Improvement over that of Ref. 5 
as originally forecast, at least for predicting these coaxial 
jet data. 

The mixing rates predicted with the TKE model and computa- 
tional technique designated by Spalding as BRASS are much too 
fast resulting in very rapid velocity decays The relative 
Influence of the walls whicn could not be moved out beyond abort 
3.5 In. without computational Instability resultlno, and the 
influence of the TKE model which overpredi cted the mixing dras- 
tically could not be assessed. 

'Adjustment" oi the velocity field, i.e., the profile 
shapes, no doubt could have been made by further modification of 
the constants of the models. However, such modification was 
beyond the scope of this work, which had as its objective an 
evaluation of existing modelr . Results clearly demonstrate that 
a new mixing model , or at least drastic revision of existing 
models. Is required before computation of practical coaxial jet 
flows in which the velocity of the outer jet Is siqnlf icantly 
greater than that of the central jet, and which contain a splitter 
plate, can be made with confidence. 
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V. CONCLUSIONS 

A simple, semi -empiri cal performance correlation/prediction 
technique applicable to gaseous and liquid propellant rocket 
engines was presented. Correlations were attained by "adjusting" 
the computation of the gaseous mixing of an unreactive, coaxial 
jet using a Correlation Factor, F, which resulted in prediction 
of the experimental combustion efficiency to within 0.1% for each 
firing. The technique was successfully applied to Rockeldyne, 
Aerojet, TRW, and Bell Aerospace gaseous H 2 /O 2 rocket engines 
utilizing coaxial, triplet, trislot, premix, and reverse flow 
injector elements, and to Bell's 6000 lb-thrust Orbital Maneuvering 
and 600 1b-th“ust Reaction Control Engines, which utilize triplet 
and unlike doublet injector elements, respectively, and liquid 
monomethyl hydrazi ne and nitric oxide propellants. 

The range of conditions over which the gaseous H 2 /O? rocket 
engines were correlated is: Injector elements 7 to 96, Chamber 
length 3.7 to 8.3 in., L* 7.7 to 57, Chamber pressure 28 to 470 
psia, "Effective" radius (R-|) 0.026 to 0.71 in., Ug ? 90 to 580 
ft/sec, Uh ? 520 to 3720 ft/sec, Vd 1.04 to 22.6, n e vn 81.6 to 
99.5%, Fi 1 m-cool ant level 0 to 30%, O 2 /H 2 mass ratio 1.9 to 7.5, 

F 0.029 to 2.2, and F * L / R i 3 to 17. For these test data, F was 
essentially independent of hexp’ num ^ er of injector elements, 
chamber pressure, and (except in one injector) film-coolant level. 

Predictions of rocket engine performance using the simple 
unreactive (cold flow), H 2 /O 2 mixing calculation require as 
input: chamber pressure (Pc), length from injector to nozzle 
tnroat (L), number of injector elements (EL), and the mass flow 
rate, injector area, and total temperature of each propellant. 

An additional input required is the Correlation Factor, F, which 
can be estimated from plots of F»L/Ri versus velocity ratio, Vr , 
for the injector configuration of interest. 

The success achieved with these simple correlations suqgests 
that key rocket design parameters are: 1) injector configuration, 

2) chamber length divided by the "effective" radius of the central 
oxidizer jet, L/R], 3) velocity ratio, Vr, (rather than the actual 
velocities themselves), and 4) prooellant type. The correlation/ 
prediction technique is useful for predicting optimum operating 
conditions for a given injector geometry, and for assessing the 
consistency of test data. Before liquid rocket and Space Shuttle 
Main Engine performance predictions can be made with confidence, 
additional liquid rocket data, covering a wider range of condi- 
tions, must be correlated. 

Static pressure, mean velocity and turbulence intensity 
ofiles in the developing region of a non-reactive (air-air) 
coaxial jet were presented in which U 0 /U-j = 2.8, typical of 
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coaxial Injector elements. Detailed data were obtained at twelve 
axial locations (extending from the nozzle exit for a distance 
of five diameters) downstream from a single element of the Bell 
Aerospace H 2 /O 2 19-element coaxial Injector. Measurements of 
mass-flux per unit area (using a constant temperature anemometer), 
total pressure, and local temperature were used In the determina- 
tion of local static pressure and velocity. 

Results showed the static pressure decreases In the vicinity 
of the nozzle exit. Although the pressure reduction was only 
0.34 psl, it substantially altered the flow. Velocities In the 
central jet near the nozzle exit, decreased initially (there was 
no central velocity core) as a result of both pressure gradients 
and viscous mixing. 

The experimental results were compared with analytical 
predictions made using eddy viscosity and turbulence kinetic 
energy mixing models and available computer codes. Comparisons 
were disappointing, the eddy viscosity model underpredicted the 
extent of the mixing, especially in the reqion within the radius 
of the central jet; whereas, Spalding's TKE model drastically 
overpredicted the extent of the mixing throughout the flow 
field. In sum, in practical coaxial injectors: 1) the effects 
of pressure gradients can significantly influence the mixing and 
combustion, and 2) present turbulent mixing models and computa- 
tional techniques are not adequate for predicting the flow in 
the developing region. 
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•Bell IRAD Model 8616 - High Pressure Gaseous H-/0- and Bell IRAQ Model 8616 - Low Pressure Gaseous H-/0 7 (Ref. 13) 
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FIG JR E 2 . Schematic of Coaxial Turbulent Jet and Definition of Mixing Regions 
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FIGURE 3C . FINAL VALUE OF TURBULENT (EDDY) VISCOSITY, e . ROCKETDYNE COAXIAL INJECTOR 































FIGURE 4B . CORRELATION OF ROCKE1 
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Figure 6 A .Aerojet Premixed Injector Performance Data 
Reference 1 0 






FIGURE 66 . CORRELATION OF AEROJET PREMIX INJECTOR PERFORMANCE DATA 
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FIGURE 6C . FINAL VALUE OF TURBULENT (EDDY) VISCOSITY, e_. ROCKETDYNE COAXIAL INJECTOR 
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Figure 7 A Beil Reverse Flow ;r.:”::tor Performance Data 
Reference 1 1 
































FIGURE 8B. CORRELATION OF TRW TRIPLET INJECTOR PERFORMANCE DATA 
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FIGURE 9A . Bell Coaxial Injector Performance Data 

Reference 1 3 
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VELOCITY RATIO V R (U H /U 0 ) 
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Figure 16. Bell Hydrogen/Oxygen Injector Element 







Figure 17. Sketch of Flow Around Splitter Plate for High Velocity 
Outer and Low Velocity Inner Jets 
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Figu.i 18. Comparison of Initial Velocity Profiles with 1/7 and 1/10 Power Laws 
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Figure 23. Centerline Velocity versus Axial Displacement 
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Figure 25. Variation of Mass Entrainment with Axial Displacement 





Figure 27. Centerline Turbulence Intensity versus Axial Displacement 
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FIGURE 30. COMPARISON OF EXPERIMENTAL AND PREDICTED MEAN AXIAL VELOCITY PROFILES, z « 0.010 In 
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FIGURE 31. COMPARISON OF EXPERIMENTAL AND PREDICTED MEAN AXIAL VELOCITY PROFILES 
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FIGURE 36. COMPARISON OF EXPERIMENTAL AND PREDICTED MEAN AXIAL VELOCITY PROFILES 



FIGURE 37. COMPARISON OF EXPERIME 





-0.60 -0.*I0 -0.20 H 

RR0IU5 (I 

FIGURE 38. COMPARISON OF EXPERIMENT 






.500 IN. 
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APPENDIX A 



GASEOUS HYDROGEN/OXYGEN ROCKET PERFORMANCE 
CORRELATION AND PREDICTION PROGRAM 


1 . Description of Program 

The computer program used for the Gaseous H 2 /O 2 Rocket 
Performance Correlation Prediction Calculations is similar to 
that described in detail in Ref. 8, i.e., it utilizes an explicit 
finite-difference numerical integration technique in which the 
governing shear layer equations (which are identical in form to 
the boundary layer equations) are solved in the von Mises co- 
ordinate system. The basic partial differential equations solved 
are: 


Mass * 


, N a N 
3puy . 3pvy 

3x 3y 


Momentum 


9u , 9u 

pu Sx + pv aF 


1 3 r N 3u-i d£ 

N 3y Lyy 3y J ~ dx 

y 


Energy 


2 2 
a u_ 3 L 

3H . 3H _ I 3 r N t_2 . 1 ,3H _ 2 ) 

pu 3x pv 3y “ N 3y ^ yy * 3y Pr 3y 3y ' 




y-1 


J 3 a 

ay 


j 


->] 


Diffusion 


(a) Species (frozen and finite rate chemistry) 


pu 


3 or 
3x 


J 


+ DV lai = 1_ L. [H_ V N isi] + ^ 

pv N 3y L Sc y 3y J 


ay 


♦Nomenclature is consistent with that of Ref. 8. 
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(b) Elements (equilibrium chemistry) 


35 k , 
pu 53T + 


pv 


3oH 

ay 


1_ l_ rU N 3a k -| 

N 9y L Sc y 9y J 

J 


v/he re 


N = 


0 

1 


Planar Two-Dimensional Flow 
Axisymmetric Flow 


The equations are transformed from the physical plane to a 
str >amfunction plane specified by 


N M M M 

9i{> = puy 9y and \p dip = pvy 9x 

In this plane, the Conservation Equations become: 

Momentum 

9u _ 1 9 r 9u -| 1 dp 

9x” - , N 9 U> *" ~ ^ v 

Energy 

8 u£ 

9H _ 1 9 r , 2 1 / 9H 

9x - dip La 1 9\|» Pr ^ 9^ 9^ ' T v Sc Pr' J £ 1 “ dip 

Diffusion 

(a) Species 

9o£ = 1_ 9_ r a_ 9ot J , w£_ 

9x 9iji L Sc 9^ J pu 

(b) Elements 

9a k _ 1_ 9_ r a_ 9a k - 1 
9x “ dip L Sc 3^ J 


0 ^ 


pu UA 


9 

3 2 N 


,1 


n 

r 


J 9a! 


1 1 


where 


2 


N 



2. Transformed Streamfunctlon Plane 


The streamfunctlon transformation may be written as 


*b u f y b N 

difi = I puy dy 

a y a 


0 - Planar 2-D Flow 

1 - Ax 1 symmetric Flow 


If the initial profile consists of two uniform streams; say 
hydrogen from yi = 0 to y£ s yn? and oxygen from y£ to y 3 , 
then the streamfunctlon can be computed as 


. N+l . N+l 

*2 - ♦, 




y? +1 ) 


But 


f] = yi s 0, so that 


For the oxygen stream 


1 

N+l 


*2 * C(pu) H yj* 1 ] 


1 

N+l 


i|» 2 = y 2^ pu ^H, 


* 3 • C*S +1 * (P«*) 0 s ,(yg +1 

1 


1 

N+l 


y! +1 )] 


Note that for planar 2-D flow, jjry a 1, the streamfunctlon is 
Identical to the mass flow rate" 


(cu) y 2 0> • *, 
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The proper way to model the Initial potential core region 
Is to put the bulk of the grid points In the stream for which 
tpuo s *P 2 or 4 * 0 ? * 4»3 - t|<2 Is smallest. Also, it should be noted 
that tp i may be^> 0 and ^+ op < 413 at the origin, and it is hence 
possible to concentrate initially on the Interface region between 
the two streams. 


The program computational grid will automatically spread 
out towards ip ^ = 0 and 1^3 as the mixing process occurs downstream 
Also, more than two streams may be modeled, as long as sufficient 
grid points are used; each stream should nave a minimum of ten 
grid points. The interval between grid points, Ai|>, at ifi = 0 is 

determined by Aip - g es -j r ed Number of Grid Points' 


3 . Brief Description of Subroutines 

The Flow Chart in Block Diagram form of the Bell Aerospace 
Combustor Correlation/Prediction Program is presented in Fig. A-l 
The important Subroutines are described briefly below: 


MAIN 

MAIN is used to input the data for each case, initialize the 
necessary variables, and then to call the true Master subroutine 
of the orogram, MARCH. MAIN will also control the FCORE varying 
of the viscosity level and repetition of a calculation until it 
converges upon experimental results. 

BEGIN 

BEGIN is used to store the fits for the thermodynamic data for 
the chemical species. 

PIF 1 

PIF 1 is used to make linear interpolations between two one-dimen 
sional arrays. 

HEAT 

HEAT is used to (a) compute the enthaloy from the specie mass 
fractions, thermodynamic fits, and temperature, or (b) invert 
to obtain the temperature from the enthalpy and specie mass 
fractions . 

DENSE 

DENSE is used to compute the mixture density, mass flux, and 
molecular weight at each point in the flow field. 



SED 

SEO contains the viscosity models and employs the Inout specified 
model . 

MARCH 

MARCH Is the master controlling subroutine of the analysis, and 
specifies the sequence of operations performed for each step In 
the axial direction. 

STEP 

STEP contains the axial step size stability criteria, and computes 
the axial step size accordingly. 

PRINT 

PRINT prints out the radial profiles of the independent and 
dependent variables at all specified axial locations. 

CONSRV 

CONSRV solves the Conservatl on .Equations for the change In the 
Independent variables (H, u, aJ or a*) across each diffusion 
step. 

GRID 

GRID adds mesh grid points, and doubles the interval between qrld 
points, as the calculation proceeds downstream. 

BULK 

BULK computes bulk properties across the flow field at each axial 
station 


4. Input Format and Deck Organization for IBM 7090 

The Input Format required for the Data Deck is shown In 
Table A-l. As many sets of data (cases) may be run as desired. 
Only one END card is needed following the last case. The 
organization of the complete deck for the IBM 7090 is Illustrated 
in Fig. A-2. A sample Input for the Bell Coaxial Injector Case 
892, ready for key punching, Is presented In Table A-2. 

5. Output Format and Sample Case 

The Output Format Is given In Table A-3 and the computer 
printout obtained for the sample case above Is shown In Table 
A-4. In this printout, ETA PRED computed at the Initial axial 
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station ( x * 0 ) and final axial station (x=l) are each given along 
with F, and ETA for each Iteration. When the value of F consistent 
with the experimental efficiency (ETA) is being determined, more 
than one Iteration Is generally required. When only a prediction 
is desired, the value of ETA must be set to zero and no Iterations 
will be made. In such cases only one computation, which terminates 
at the axial station just beyond x*L will be made, and the values 
of ETA PRED are consistent with the Input value of F. Note, that 
when Iterations are required to determine the appropriate value 
of F, the number of the iteration printed out will not be changed 
in cases in which values of ETA PRED are greater than 99.99%. 

This procedure simplified considerably the counting procedure and 
rarely occurs. 


6 . Computations and Results 


The calculations were always begun assuming slug (st^p) 
profiles and Scj = 0.7, Pr T = 1.0. The initial Yq« profile used 
to start the computation (x=0) for Rocketdyne CaseH2H, Ref. 9, 
is plotted in Fig. A-3 along with the Yq ? profiles computed at 
intermediate and final stations of 0.5, ^2.5, and 5 in., which 
were computed using the unreactive mixing calculation discussed 
in ANALYSIS. Note, that for the binary Hp/Og system considered, 

Y^p ? 1 - Yq?* so that the corresponding Yuj values can be obtained 
from the plots as well. The velocity profiles at these same 
stations are plotted in Fig. A-4. The initial step velocity 
profiles at x=0 are simply the bulk mean velocities computed using 
Eq. (9). The Yh 2 q profiles computed using either Eq. (5) or (6) 
are plotted in Fig. A-5 at the same axial stations. As 
explained in ANALYSIS the water was merely computed from either 
Hp or O 2 concentrations (whichever was the limiting reagent), and 
played n£ role in the mixing calculation; of course, the Yh« 0 
profiles were required for the computation of the combustion 
efficiency, np re d» using Eq. (8). Although more sophisticated 
computations , e. g . , including the finite rate chemistry, are 
within the capability of the computer program, the success 
achieved with the very simple approach suggested that expenditure 
of additional computer time was unwarranted - as long as only the 
combustion efficiency was to be predicted rather than details of 
the flow. 

The H 2 /O 2 mixing computations are well illustrated in Figs. 

A-3 to A-5. There is gradual mixing of both mass and momentum 
(stagnation temperature remains constant). As mass mixing pro- 
ceeds, the computed water level grows, and hence, npred con- 
tinually increases [Eq. (8)]. Startinq with a given set of ini- 
tial conditions, there is a one-to-one correspondence between the 
Yh 2 0 profiles and npred* In the iteration to determine the ap- 
propriate value of F, the profiles at the final axial station, 
x»L, are "adjusted" until npred computed from Eq. (8) yields nexp 
to within the desired tolerance, e.g., 0.1%. 
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Values of the turbulent viscosity, l.e., eddy viscosity, 
e n , used throughout the mixing computation for Rocketdyne Tests 
12H, 14H , and 18H, Ref. 9, are plotted In Fig. A-6. Note that 
by assumption e n Is not a function of radial position, and so 
Is constant at each axial position. These particular test cases 
were selected because they covered the entire Vp range of the 
Rocketdyne Coaxial Injector test data; Test 14H is at the maximum 
of the curve of F vs Vo In Fig. A-7. For comparison, values of 
enjodel co,n P ute<1 directly from Eq. (4) are presented In Fig. A-8. 
Of course, these values are directly proportional to the e n 's 
plotted In Fig. A-6 [Eq. (1)]; however, their relative magnitudes 
are shifted because the F's vary (as do the n's). In view of the 
erratic behavior of these mixing coefficients. It is very sur- 
prising that useful correlations, such as those presented In 
Fig. A-7 were obtained; however, useful correlations have been 
obtained for every case to which the technique has been applied. 
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Table A-l 

Bell Aerosoace Combustor Correlation InDUt Format 
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Length From Injector to Nozzle Throat 
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Table A-3 

Bell Combustor Correlation and Prediction Program 

OUTPUT FORMAT 




LINE 

NAME 

FORMAT 

1 

PROGRAM 



HEADING 


2 

TITLE 

18A4 

3 

ETA 

F15.6 


F 

F15.6 


PC 

F15.6 


L 

F15.6 


EL 

F15.6 


VRINJ 

F15.6 


O p !NJ 

F15.6 


R1IN 

F15.6 

4 

MD0TH2 

F15.6 


AH2INJ 

F15.6 


TH2 

F15.6 


MD0TH2/EL 

F15.6 


AH2INJ/EL 

F15.6 


UH2INJ 

F15.6 

5 

MD0T02 

F15.6 


A02INJ 

F15.6 


T02 

F15.6 


MD0T02/EL 

F15.6 


A02INJ/EL 

F15.6 


U02INJ 

F15.6 

6 

X 

F15.4 


L 

F15.2 


F 

F15.6 


ETA PRED 

F15.3 


ETA 

F 1 5. 1 


ITER 

19 


DESCRIPTION 


Input Run Designation, ID, etc. 

Experimental Combustion Efficiency (%)* 
Input Correlation Factor 
Input Chamber Pressure (psla) 

Input Chamber Length (In.) 

Input No. of Injector Elements 
Velocity at H2/Veloc'1ty 02 Injector 
O 2 /H 2 Mass Flow Ratio at Injector 
Radius of Central Jet (effective radius 
If not axlsymmetrlc Injector) (In. ) 


Input l ifi H at Injector (lbm/sec) 

2 2 

Input l Area H 2 at Injector (In. ) 

T H 2 , total <° R) 

7 ifi u /Injector Element (lbm/sec-element) 

H 2 2 

[ Area H 2 /Injector Element (In. /element) 

H 2 Jet Velocity (in annulus ) (ft/sec) 

Input l ifi Q Injector (lbm/sec) 

2 2 

Input l Area 0 2 Injector (In. ) 

Input T o 2 ,total ^ 

l /Injector Element (lbm/sec-element) 

2 2 

1 Area 0 2 /lnjector Element (In. /element) 

0 2 Jet Velocity (central jet)(ft/sec) 


Axial Station ( In . ) 

Input Chamber Lenqth (In.) 

Input Correlation Factor 

Computed/Predicted Combustion Efficiency {%) 
Experimental Combustion Efficiency (%) 
(repeated for comparison) 

No. of iteration when Input ETA>0; other- 
wise blank 


♦If ETA-0 predict Combustion Efficiency using Input value of the 
Correlation Factor, F. If ETA>0, e.g., 95%, a value of F con- 
sistent with the input ETA will be determined by Iteration. 
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Figure A-3. Mass Fraction Oxygen at Various Axial Stations Computed using of Figure A-6 

Rocketdyne Coaxial Injector Test 1 2H, Ref. 9 
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Figure A-4. Velocity Profiles at Various Axial Stations Computed using e„ of Fig. A-6 
Rocketdyne Coaxial Injector Test 1 2H, Ref. 9 
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Figure A-5. Water Profiles at Various Axial Stations Computed using Eqs (5) and (6) - 

Rocketdyne Coaxial Injector Test 12H, Ref 9 
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TURBULENT VISCOSITY, 



AXIAL DISTANCE, z. IN. 

Figure A-6. Value of Turbulent (Eddy) Viscosity used in rj Computation for Rocketdyne 

Coaxial Injector Tests 12H, 14H, and 18H. Ref. 9 
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Figure A-8. Comparison of Predictions of Mixing Model, [Eq. (4)] for Rocketdyne 
Coaxial Injector Tests 12H, 14H, and 18H, Ref. 9 
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NOMENCLATURE 


A Area, ft 2 

2 

£A Total Injection area for a particular propellant, ft 

. .1 5 

D p Empirical density parameter (= P u /P m j n ) 

EL Number of injector elements 

F Correlation Factor (= e /e n)0(Jel ) 

f Empirical_length parameter (h 0.5 + 0.005F; F < 100, 
and 1.0; z > 100) 

G Total mass flow rate of gaseous propellant measured at 
injector, lbm/sec 

L Length of combustion chamber (injector to nozzle throat), 
ft 

M Mach number 

MW Molecular weight, lbm/lb-mole 

rti Comouted mass flow rate, lbm/sec 

2 

P Static pressure, lbf/ft 

P_ Chamber pressure, lbf/ft 2 

w 

2 

?' Static pressure during hot film calibration, lbf/ft 

w 

2 

P Q Total pressure, lbf/ft 

Pry Turbulent Prandtl number 

R Universal gas constant, 1545 f t- 1 bf /I b-mol e»°R 

R, "Effective" radius of central oxidizer Injector, 

1 defined in Eq. (11), ft 

R, "Effective" radius of annular fuel Injector, 

* defined in Eq. (12), ft 

r Radial coordinate, ft 

Scy Turbulent Schmidt number 

T Static temperature, °R 
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T 0 Total or stagnation temperature, °R 

U Axial velocity, ft/sec 

U Time mean (average velocity, ft/sec 

/IF Turbulence Intensity, ft/sec 

V Indicated mean velocity (hot film), ft/sec 

V R Velocity ratio (= U H /Ug ) 

Y Mass fraction 

z Axial coordinate, ft 

z Normalized axial coordinate (= z / 2R ^ ) 
y Ratio of specific heats 

e model Turbulent (eddy) viscosity predicted by Eq. (4), 

1 bm/sec-f t 

e Turbulent (eddy) viscosity used In computation of 

Ved (i Fe model > ’ ' bm/sec-f t 

n avn Experimentally determined combustion efficiency 

V A P 

n pred Combustion efficiency predicted using computer program 

3 

p Density, Ibm/ft 

3 

p. Density during calibration of hot film, Ibm/ft 

C 

Subscripts 

e External or freestream conditions 

1 Central or Inner 

L Splitter plate or land 

o Annular or outer 

u Condition at velocity half radius where U=0.5(U m „„ + 

U ) max 

min' 

» Freestream condition 
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